Charge Order Superstructure with Integer Iron Valence in Fe20B03 
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Solution-grown single crystals of Fe20B03 were characterized by specific heat, Mossbauer spec- 
troscopy, and x-ray diff'raction. A peak in the specific heat at 340 K indicates the onset of charge 
order. Evidence for a doubling of the unit cell at low temperature is presented. Combining struc- 
tural refinement of diffraction data and Mossbauer spectra, domains with diagonal charge order are 
established. Bond-valence-sum analysis indicates integer valence states of the Fe ions in the charge 
ordered phase, suggesting Fe20B03 is the clearest example of ionic charge order so far. 

PACS numbers: 61.50.Ks, 71.30.-fh, 71.28.-|-d, 61.10.Nz 
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Many physical phenomena in transition metal ox- 
ides, including colossal magnetoresistance [H and high- 
temperature superconductivity 0], are related to charge 
ordering (CO). Ideally, CO consists of charge carriers lo- 
calizing on ions with different integer valences forming 
an ordered pattern However, the application of this 
"ionic CO" concept has been controversial 0,01, because 
observed valence separations are usually small, and there 
is no clear example of CO with integer valences. 

Mossbauer spectra on the mixed-valent warwickite 
Fe2 0B03 suggested a large, though not quantified, Fe 
valence separation below the onset of a nionoclinic dis- 
tortion of the structure (Fig. [l]) at 317 K j6|. It is natu- 
ral then to suspect an ordered arrangement of Fe^"*" and 
Fe'^^ ions (ionic CO), and Fe20B03 has been suggested 
as an example of electrostatically driven CO 0. How- 
ever, no experimental evidence of a CO superstructure 
was found on the available polycrystalline samples, and 




FIG. 1: (Color online) Crystal structure of Fe20B03 at 
355 K. The unit cell (shaded yellow) is orthorhombic (Pmcn), 
at lower T the symmetry is lowered to monoclinic, with a 
CO superstructure having a 2axbxc cell. Structurally dis- 
tinct FelOe and Fe206 octahedra build ribbons of four edge- 
sharing chains (numbered, see text) along a. Detail views of 
the ribbons at 355 and 100 K are given in Fig. |4] 



consequently the occurrence of CO in Fe2 0B03 has been 
under debate 

Here, we report the first observation of superstructure 
reflections in single-crystalline Fe2 0B03, using X-ray 
diffraction. Combining structural refinement, Mossbauer 
spectroscopy, and electronic structure calculations, we es- 
tablish a diagonal CO configuration. Bond-valence-sum 
analysis indicates that the ordered iron valence states 
are very close to integer Fe^+ and Fe'^^. Thus, Fe2 0B03 
is an excellent example of ionic CO. We discuss implica- 
tions of the large structural modulations on the relevance 
of the electron-lattice coupling in driving the CO. 

Needle-like single crystals of Fe2 0B03 (Fig. [2] inset) 
with length up to 1.5 cm were grown from a flux with a 
procedure very similar to the growth of Fei.giVg.ogOBOs 
reported by Balaev et al. [8|, except that we omitted 
V2O3 from the flux to avoid V doping. ^^Fe Mossbauer 
spectra, obtained on powdered crystals using a constant- 
acceleration spectrometer Q, were similar to previous 
results [l3| with isomer shifts at low T (Fig. [Sb) indi- 
cating divalent and trivalent Fe with no electron hopping. 
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FIG. 2: (Color online) Specific heat Cp/T, measured by ac 



calorimetry. Lower inset: 
Fe20B03 crystal. 



around 340 K. Upper inset: 
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TABLE I: Refinement results [IJI for Fe sites in Fe20B03 at 
355 K (space group Pmcn, a = 3.18 A, & = 9.40A, c = 9.25A) 
and 100 K {average structure, space group Pc, 2a = 6. 33 A, 
b = 9.38 A, c = 9.25 A). Average Fe-O bondlengths, bond- 
valence-sum, and distortion parameter (see text). 



Site 


{d(Fe-0)>/A 


Bond Valence Sum 


Distortion F 


355 K 


Fel 


2.094(2) 


2.39(1) 


2.8 


Fe2 


2.088(2) 


2.46(1) 


4.0 


100 K 


Fela 


2.140(6) 


2.04(3) 


5.9 


Fe2a 


2.154(6) 


2.03(3) 


7.5 


Felfc 


2.020(6) 


3.02(5) 


0.4 


Fe2b 


2.019(6) 


3.07(5) 


0.0 


Felc 


2.098(7) 


2.35(8) 


2.9 


Fe2c 


2.087(6) 


2.47(8) 


3.2 


Feld 


2.098(6) 


2.34(8) 


2.9 


Fe2d 


2.086(6) 


2.48(8) 


3.6 




FIG. 3; (Color online) Composite X-ray diffraction precession 
images with intensities extracted from about 500 individual 
frames, a: 355 K. Due to an n glide plane, spots with h + k 
odd are systematically absent, b: 100 K. Weak additional 
spots (one indicated by an arrow) index to {h+^,k,0), indi- 
cating a superstructure with doubled (2a x b x c) unit cell. 



The specific heat (Fig.[2|) was measured with ac calorime- 
try ll| on several crystals giving consistent results. A 
step- like feature around 153 K is due to the known mag- 
netic transition [l^]. A peak at 340 K with '^0.25K hys- 
teresis suggests an additional, weakly first-order, phase 
transition, which according to a powder diffraction study 
[ist corresponds to the monoclinic-orthorhombic transi- 
tion, and which we attribute to to the onset of CO. The 
estimated entropy change associated with the transition 
is very small (~0.01i?/Fe ion, R = 8.14 J/mol/K), but 
we point out that ac calorimetry is sensitive only to the 
reversible heat flow. 

To corroborate the CO, we collected single-crystal X- 
ray diffraction data with a Bruker SMART APEX CCD 
diffractometer, using MoKa radiation. Figure [3] shows 
composite precession images of the hkO plane with inten- 
sities extracted from about 500 individual frames mea- 
sured at 355 K (a) and 100 K (b). The pattern at 355 K 



is expected for the Pmcn (no. 62) space group, and the 
refined structure (Figs. [1] and fl4j) was very similar 
to the one reported for poly crystalline Fe20B03 0. At 
100 K, weak additional spots indexing to {h+^,k,0) in- 
dicate a superstructure with 2axbxc cell, attributed 
to CO as detailed below. Visual inspection of the raw 
data indicated overlapping peaks of roughly equal in- 
tensity consistent with monoclinic twinning with similar 
weight of domains with opposite sense of monoclinic dis- 
tortion. The twinned peaks were not sufhciently resolved 
and therefore the refinement was conducted for the com- 
bined peaks. As a result, the monoclinic angle f3 could 
not be refined, and was set to 90°. 

To assess possible symmetries, the data were first re- 
fined with no symmetry constraints (space group PI). 
A Pc (No. 7) space group symmetry was apparent, 
and a corresponding refinement (Fig. had a resid- 
ual R[F'^ > 4ct(F2)] = 4.62% significantly better than 
alternative space groups T4i]. Small (very weak inten- 
sity) reflection condition violations and a Flack parame- 
ter close to 50% suggest that the structure is not homoge- 
neous, but an average of domains with different, possibly 
centro-symmetric, structure. Room temperature diffrac- 
tion data are very similar to the 100 K data, indicating 
negligible influence of the magnetic transition on the CO. 

For octahedrally coordinated Fe^+ and Fe'^"'" ions, 
the expected average Fe-0 bondlengths are 2.16 and 
2.02 A, respectively 15 [. For the refined structure, Fe-0 
bondlengths for half of the eight Fe sites are intermediate 
and similar as at 355 K (Table suggesting CO in only 
half of the chains. However, Mossbauer spectra show that 
the local structure is different because there is no interme- 
diate valence Fe in the low T phase (Fig.[S)3). The refined 
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355 K 100 K ; average structure = "up" - domain + "down" -domain 




FIG. 4: (Color online) Charge order in the four-chain ribbons (numbered as in Fig. [T}. a: 355 K. Atoms are drawn as thermal 
ellipsoids. For one FeOe octahedron a local coordinate system is indicated. b,c: CO at 100 K (high/low valence Fe is shaded 
dark/bright): The "global" structure, refined in Pc (b) arises from averaging domains with different diagonal order (c). Large 
oxygen shifts within the superstructure are emphasized by red arrows in b. 



structure thus arises from averaging local structures of 
two types of domains. In two of the chains, II and IV in 
Fig. [31 structural distortions are averaged out, indicating 
that these sites have opposite valence, and thus distor- 
tions, in the two domain types. The CO in the other 
two chains is clearly preserved, indicating that their Fe 
sites have identical valences in all domains. Because the 
Fe sites with equal a position in the chains with globally 
preserved CO have opposite valence, the local CO con- 
figuration is diagonal, with "up" and "down" diagonals 
making up the two types of domains (Fig.|3j;). 

To corroborate the local diagonal CO, we per- 
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FIG. 5: (Color online) Integer valence separation, a: Fe va- 
lence from bond-valence-sum analysis on four sites indicated 
in Fig. |4]d (-I-), compared to the classical CO example mag- 
netite (B sites, o[iy|) and the clearer example EuBaFe205 
(■[l3|). b: Mossbauer isomer shifts relative to aFe (-I-), com- 
pared to magnetite (B sites, o[l^) and EuBaFe205 (■[l3|)- 
The two horizontal lines represent isomer shifts [T?] for Fe^+ 
in FeO (blue) and Fe^^ in Fe203 (red) at room temperature. 



formed first principles electronic structure calculations 
(GGA-I-U, with fixed cell, but in contrast to 20] opti- 
mized atom positions), using techniques as described in 
'2^. The GGA-fU calculations (C/ = 5.5eV, J = 0.89eV 
20|) started with the crystal structure as reported by 
Attfield et al. ^ with doubled a axis, in which all Fe 
atoms have a uniform valence Fe^'^+. As in (iSl, a diag- 
onal configuration with charge separation was obtained 
even with fixed atom positions. Relaxing atom positions 
led to an additional decrease in energy (total CO gain 
174.5 meV/Fe ion), with CO distortions of the optimized 
atom positions (shown in Fig. Uj;) qualitatively equal to 
the refined ones (neglecting the averaged out sites in 
chains II and IV), though about 20% less in magnitude. 

The shifts in the oxygen positions (Fig. arrows) 
cause the average Fe-0 bondlengths for sites Fela and 
Fe2a to increase and for sites Felb and Fe26 to decrease. 
Valence states and ion sizes are intimately related, allow- 
ing the length of the Fe-0 bonds to be used to calculate 
the Fe valence through bond-valence-sum (BVS) analysis 
fH. The BVS is 



y = ^ exp [(do - di) /0.37] , 



(1) 



where V is the valence of an ion to be determined, di 
are the bondlengths to other ions, and do is a tabulated 
[2^ 1 empirical parameter characteristic for a cation-anion 
pair. At all temperatures the BVS of B and O are close to 
3 and —2, respectively, as expected. A complication for 
Fe is that the empirical do are slightly different for differ- 
ent valence states. Using tabulated [2^] do for both Fe^"*"- 
O and Fe^+-0, the calculated valences for ions Fela,Fe2a 
and Fel6,Fe26 are within 0.2 of 2 and 3, respectively. As 
is customary in this case [23| , the final valences were then 
calculated using do of Fe^^-O for ions Fela and Fe2a and 
do of Fe'^"'"-0 for ions Fel6 and Fe26. For the averaged 
sites in chains II, IV at low T and all Fe sites at high T the 
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BVS are intermediate between 2 and 3, and the final va- 
lences were obtained by averaging the valences calculated 
with do for Fe^+ and Fe^+. 

Averaging over anti-phase domains can only decrease 
the difference between high and low valence Fe-0 
bondlengths in the average structure. Consequently the 
valence separation from BVS obtained from any refined 
global average structure provides a lower limit for the ac- 
tual local separation. We focus on the sites in chains I 
and III, which provide the most stringent limit. These Fe 
valences are integer within the resolution of the method 
(Fig. [Hi). FezOBOa is the first CO oxide for which BVS 
clearly indicate integer valence for both valence states. 
The valence separation is considerably larger than in 
the so far clearest examples, YBaFe205 ji^l and related 
compounds (^^0.7), and much larger than in the classi- 
cal, though not fully understood, CO example magnetite 
(<0.4 [16|) or in the colossal magnetoresistance mangan- 
ites (~ 0.45 [i]). The large valence separation suggests 
that Fe2 OBO3 is an ideal example of ionic charge order. 

In agreement with the large valence difference, 
Mossbauer spectral isomer shifts S (Fig.O)) for two dou- 
blets are close to i5 of the Fe'^+ compound Fe203 and 
for the other two doublets to S of the Fe^+ compound 
FeO. The reason that Fe20B03 is unique among oxides 
in demonstrating ionic CO may be attributed to an an- 
tagonistic inductive effect |25l| from the B-0 bonds. Since 
B is more electronegative than Fe, O prefers to share its 
electrons with B, making the Fe-0 bonds more ionic. 

The large (> 0.2 A) oxygen shifts along a accompany- 
ing the CO also affect the distortions of the coordination 
octahedra. At high T the largest distortion is an elonga- 
tion along a (Fig. 2^), leading to a distortion parameter 
r, defined as the difference in % between bondlengths 
in the local xy plane (Fig. [4^) and those perpendicular 
to it, of 3 to 4%. In the CO state, this distortion is re- 
duced almost completely for Fe^+ sites, but increased for 
Fe^^ sites, consistent with the extra electron occupying 
a dxy orbital with lobes \\a. Because the distortion lifts 
the degeneracy of the t2g orbitals, an energy gain associ- 
ated with the Jahn- Teller effect |26j likely contributes to 
the overall energy gain of charge localization and order. 
Thus, CO in Fe20B03 is not driven entirely by electro- 
static interactions between carriers, although the latter 
alone seems to be sufficient in establishing CO, as indi- 
cated by the GGA-I-U calculations with atom positions 
fixed. This exemplifies the relevance of electron-lattice 
effects in correlated oxides, which likely is ubiquitous. 

In summary, superstructure reflections at in Fe2 0B03 
arise from diagonal CO with two domains and Fe valence 
states very close to integer, suggesting that Fe20B03 is 
an ideal example for ionic CO. Magnitude and form of the 
atomic displacements in the CO suggest that electrostatic 
energy is not the only relevant energy scale, but coupling 
to the lattice is important as well. Very anisotropic dis- 
placement parameters of the oxygen atoms already at 



355 K (Fig. |4^) hint at significant precursor effects to the 
CO transition; a detailed study of the evolution of the 
CO with temperature is in progress. 
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